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Abstract 

The possibility that the epsilon expansion can predict the order of phase transi¬ 
tions in three dimensional field theories is examined. For a Hermitean matrix-valued 
order parameter, the epsilon expansion predicts fluctuation induced first order phase 
transitions. We analyze two 2+1-dimensional quantum field theories which exhibit 
spontaneous symmetry breaking and have martix order parameters. Using the large 
N expansion, we show that these models exhibit second order transitions and discuss 
the implications for the chiral symmetry breaking transition in 2+1-dimensional QCD 
for a critical number of quark flavors. 


1 Introduction 


In many of the applications of quantum field theory where phase transitions play a role, 
the order of the transition is often an important question. For example, the generation of 
nonzero baryon number in the early universe requires, besides interactions violating B and 
T, that the system is out of equilibrium during the process. This happens, most naturally, in 
conjunction with supercooling and reheating phenomena near a first order phase transition. 
On the contrary, a second order chiral phase transition is needed in theories of dynamical 
mass generation, and also in composite Higgs theories, to ensure that the vast array of 
mass scales generated dynamically lie well below the scale of interactions responsible for 
the relevant dynamics [II]. The nature of the chiral phase transition in QCD will also affect 
results of impending heavy ion experiments. 

There has been a substantial amount of work on understanding the detailed nature of the 
chiral symmetry breaking phase transition both in four-dimensional and three-dimensional 
QCD. Though there are some important differences, the latter, more recently studied case 
of three-dimensional QCD shares many of the features of its physical, four dimensional 
relative. The question of dynamical mass generation can be framed in much the same 
way as in four dimensions and the technical complexity of the problem is very similar. 
Many of the analytical methods which are used in four dimensions can also be used in 
three dimensions. This gives a nontrivial test of the methods in a slightly different context. 
They can furthermore be compared to numerical data which, though still very incomplete 
(particularly with dynamical fermions), should be easier to obtain in a lower dimension. 

Lower dimensional gauge theory can also be of physical interest, as in phenomenological 
models of certain condensed matter and statistical systems. The existence or non-existence 
of, and the detailed nature of phase transitions is of great importance for those models. 

A variety of methods have been used to investigate the nature of dynamical mass genera¬ 
tion in 3 dimensional gauge theories. One of the most important is the solution of Schwinger- 
Dyson equations in the quenched ladder approximation. Another very important method 
is the e-expansion of Wilson and Fisher |2j around d — 4. (We shall call this the “4 — e 
expansion.) Other methods include large N expansions, strong coupling expansions, and 
the e expansion about two dimensions applied to chiral Lagrangian sigma models (which we 
shall call the 2 + e-expansion). 
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In this paper we would like to discuss some of the relationships between these methods. 
Indeed, in the current literature there is an apparent contradiction between the analysis of 
the 4 — e-expansion and the use of Schwinger-Dyson equations. We shall comment on some 
of the issues involved. 

2 Three-dimensional QCD 

Consider the model defined by the action 

S = j d 3 x ^7 • Dip + ^j^trF^ , ( 1 ) 

with N c colors and Nf flavors of quarks: a — l...N c ,a — 1... Nf. The Dirac matrices 

are 2x2 and Hermitean and D = d + iA , F /IU = d^A v — d v A^ + i[An, A v \. 

There is no chiral symmetry in 3 dimensions. Nevertheless, mass terms can violate parity 
and time reversal invariance. Parity transformation is defined as ip(x) —> a 2 'ip(x'), for the 
fermions, and A^x) —> P^A^x') for gauge bosons. Here x^ = P^x^ where Pi = —P 2 = 
P 3 = 1 . The mass term, is obviously odd under parity, but preserves U(Nf ) global 

symmetry. 

It is known that parity cannot be spontaneously broken in a vector-like gauge theory 
. For even Nf one can have a mass term that is parity invariant, but breaks the flavor 
symmetry to SU(Nf/ 2) x SU{Nf/2 ) x 17(1). Such a mass term is of the form 

Lm 1pQN f 1p, ( 2 ) 

where Q^ f is a diagonal matrix with a set of Nf/2 eigenvalues of 1 and a set of Nf/2 
eigenvalues of -1. Then we can define a parity operator which exchanges fermions pairwise 
between the two sets. 

Whether the flavor symmetry is broken is a dynamical question. The symmetry breaking 
should generate an expectation value of a mass operator and the physical fermions of the 
resulting model would be massive. This question is complicated by the fact that, as in 
four dimensional QCD, the model likely confines color charge. It has been speculated that 
if chiral symmetry is not broken, even this three dimensional theory, in spite of its severe 
infrared behavior, would not be confining - the screening of the strong long ranged attractive 
interactions which would be necessary to restore chiral symmetry would also be sufficient to 
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relax confinement ||. As yet, this heuristic picture has little support from rigorous analytical 
calculations, but is consistent with all known analytic and numerical results. 

The question of whether there can be a phase transition between chirally symmetric and 
non-symmetric phases and if there exists such a phase transition, what is its order, has 
attracted wide interest. First, we observe that this question is unconventional in the sense 
that we are not considering the model (flj) at finite temperature, as would be the case in 
ordinary critical phenomena. Instead, we are studying how the properties of the ground 
state of ([[[) change as we vary the relevant parameters, i.e. the so-called quantum critical 
phenomena. Second, we observe that, in the conventional sense, the action ([[]) does not 
contain any coupling constants which can be varied in order to produce this phase transition. 
The QCD coupling, e 2 , is dimensional and serves to cut off the ultraviolet divergences [@. 
Its presence renders the model super-renormalizable. We therefore do not expect any phase 
transition as we vary e 2 - it is simply a dimensional parameter which fixes the units by 
which we measure energy and distance. The only parameters which can be changed are the 
numbers of colors and flavors, N c and Nf respectively. We ask the question, does the ground 
state of ([[;) exhibit different phases for different N c and Nf? Furthermore, if we vary these 
continuously (which of course cannot be done in the physical world since they are integers), 
what is the order of the phase transition? 

3 Review of investigations of symmetry breaking 

In this Section, we shall review some of the approaches to the dynamical symmetry breaking 
problem for three dimensional QCD. Generally, they are based on three different limits of 
the theory: 

• the limit where N c » Nf where planar diagrams dominate 

• the limit Nf » N c where the bubble diagrams of the large Nf expansion dominate 

• the conventional lattice strong coupling limit e 2 /A —> oo where A is the ultraviolet 
cutoff 

The most thoroughly analyzed limit of this model is that with Nf » N c . In that limit, 
internal gluon exchanges are suppressed by factors of N c /Nf and the gauge interactions are 
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very similar to those in the Abelian three dimensional electrodynamics. There have been 
extensive studies of the symmetry breaking problem in the 1 /Nf expansion using the self 
consistent Schwinger-Dyson equations || [5]. In the leading order, this approximation sums 
both bubble and rainbow diagrams and obtains j3| 0 the following picture: At large Nf 
the symmetry is unbroken if Nf > Nf Tlt , where Nf Tlt = (64/37r 2 )(A r2 — l)/iV c and is broken 
if Nf < Nf rit . The phase transition at Nf Tlt is of second order. Nf rit is large enough that 
the large Nf approximation is thought to be quite accurate in that region and corrections 
to Nf Tlt , for example, have been shown to be small in the next-to-leading order. 

The limit N c » Nf has seen much less analysis. In this limit, all of the usual planar 
diagrams of the conventional large N c expansion of QCD contribute to any process. This 
limit is therefore more likely to exhibit the features of a strongly interacting gauge held 
theory - chiral symmetry breaking and confinement. 

If the latter limit of the theory is indeed chirally non-symmetric and confining, it could 
be investigated by effective lagrangian methods. The quantum fluctuations in the broken 
symmetry phase are described by the Grassmanian SU(Nf)/SU (Nf/2 ) x SU(Nf/2). It can 
be parameterized as the set of Nf x Nf hermitean matrix-valued fields (j>(x) which have Nf/2 
eigenvalues 1 and Nf/2 eigenvalues —1. Algebraically, this can be expressed as the condition 
that 4> 2 (x) = 1. The effective action is the Grassmanian sigma-model || 

S'efj = [ d 3 x—~trdcp ■ d(j) + Nc^n [ d 4 xtr0(d0) 4 (3) 

J 2 F Jb 4 

where we have added the final, topological term and where B 4 is a four dimensional space 
whose boundary is the three dimensional space on which the first term is defined. As has 
been pointed out by Feretti and Rajeev ||, the presence of the topological term is necessary 
to break an extra, unwanted symmetry of the first term (p(x) —> — 4 >(x) which is not a 
symmetry of QCD. Furthermore, it is also necessary in order to obtain the correct current 
algebra. We observe that the number of colors, N c can affect this model only through the 
dependence on the pion coupling constant F on N c and the coefficient of the topological 
term. 

This model is difficult to analyze directly. It is not renormalizable and a large Nf ex¬ 
pansion corresponds to the as yet intractable summation of planar fat graphs. However, 
the coupling constant F is naturally of order Nf/N c , so if the number of colors is much 
larger than the number of flavors (in order to suppress the planar graphs with many internal 
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0-lines), this model should be described by the weak coupling limit. It has recently been 
studied in the saddle point approximation J7|] with results suggesting the presence of a second 
order phase transition to the symmetric phase at some critical value of Nf/N c . 

Furthermore, the model without the topological term has been studied in the 2 + e 
expansion by Brezin, Hikarni and Zinn-Justin jSj. They showed that it has an ultraviolet 
stable fixed point at a non-zero value of the coupling constant F. In ref. [7|] this was 
interpreted as implying a critical value of N c /Nf above which the chiral symmetry is broken 
and below which the chiral symmetry is restored. 

These results are compatible with the general picture of symmetry breaking which we 
advocate. However, as we shall discuss later, they may not give a complete description of 
QCD in the present context in that they neglect the presence of massless fermions at the 
phase transition point. We shall return to this point later. 

The phases of three dimensional QCD have also been studied in the strong coupling limit 
of lattice gauge theories . It was found that the strongly coupled lattice gauge theory is 
equivalent to a SU(Nf/2) Heisenberg antiferromagnet with spins taking values in the Lie 
algebra []T0[ 


h afm = e 


( 4 ) 


lattice 


where x and y are nearest neighbors on a bipartite lattice, a — 1, ..., (Nf/ 2) 2 — 1 and 


S a (x),S b (y ) = if abc S c (x)8(x-y) 


( 5 ) 


is the Lie Algebra of SU(Nf /2). The spins are in the representation given by the rectangular 
Young tableau with Nf/ 4 rows and N c columns. 

Furthermore, it was argued that the Neel order parameter of the antiferromagnet is 
equivalent to the mass parameter of the lattice gauge theory. Large N c is the classical limit 
where the ground state has Neel order. The small N c /Nf limit is the quantum limit for 
the antiferromagnet and is though to be disordered, at least in dimensions less than four. 
However the precise nature of the ground state in that case is not known. Between these 
limits there is a phase transition which separates the ordered and disordered phases. This 
transition is believed to occur at a critical line Nf ~ N c and to be of second order. It 
has been conjectured that this phase transition coincides with the chiral symmetry breaking 
transition in three dimensional QCD [7|. 
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There are numerical studies of U{ 1) gauge theory with Nf flavors [[□]]. If the phase 
transition is indeed independent of the color group, these investigations are also relevant. 
They also provide a critical value for the number of flavors. 


4 Effective bosonic theory and 4 — e-expansion 


Pisarski and Wilczek |T2j] and Wilczek [|13|] used analogies with condensed matter systems 
and the 4 — e-expansion around dimension 4 of Wilson and Fisher f2| to argue that the finite 
temperature chiral phase transition of massless four dimensional QCD is second order for 
two flavors and first order for more than two flavors of quarks, irrespective of the color gauge 
group. The first order transition is fluctuation driven and it changes the simple mean field 
theory result that one obtains using the Landau-Ginsberg effective action. In the literature 
this phenomenon is called the Coleman-Weinberg instability [|14[ . 


Pisarski has further argued [|Pj that a similar situation arises in zero temperature three 
dimensional QCD. In the following we shall review his reasoning. 

First of all, it is assumed that the relevant fluctuating degrees of freedom are bilinear 
combinations of fermion fields averaged over color. They will be denoted by Xafs {' x )> defined 
as 


N c 


X*p( x ) = Y,W X )^a( X ))' 


a= 1 


The effective Landau-Ginzburg action for such bilinears in 4 — e dimensions is 


Slg = / d '- e x\\ Tr(Vy Vy) + 


87 T Z /i 


2 ,,e r 


9i(^X 2 ) 2 + ^TrX 4 


( 6 ) 


( 7 ) 


v 2 " v 4 ! 

Note that the effective theory is independent of N c . Though the trace of the right hand side 
of (|j) does not vanish boson x in (0) is defined to be traceless. Based on general arguments, 
the symmetry breaking pattern is such that parity is conserved. |TB| That is automatically 
satisfied for a boson in adjoint representation, but not for the SU(Nf) scalar. Consequently, 
the SU(Nf) scalar does not participate in the critical dynamics. 

Pisarski |l5j used 4 — e expansion to investigate the phases of the effective theory. The 
beta functions of the effective theory are 


IV?+ 7 9 21V? 

+ —91 + 2 f 


3 N 


f 


5 N 2 f + 3 , 
-9i92 + -^2-92 


P 2 = ~eg 2 + 2gig 2 + 2 - 


IV?-9 


3 N 


f 


-9 2 - 


( 8 ) 
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There is an ultraviolet stable fixed point at the origin. There are infrared stable fixed 
points in the perturbative regime only when Nf = 2. (When Nf < \/5 ~ 2.3 there is 
an infrared stable fixed point at g* = ( 6 e/(iV| + 7),0) where the model (^) is in the same 
universality class as the 0(Nj — 1) non-linear sigma model.) When Nf > 3 there are no 
infrared stable fixed points. The situation is very similar to the one encountered by Coleman 
and Weinberg jnj]. There, the existence of two competing coupling constants allowed the 
existence of a minimum of the effective potential in the perturbative regime, leading to a 
fluctuation induced first order phase transition. Conditions for the existence of first order 
phase transitions can be obtained from the /3-functions alone as it has been shown a long 


time ago [[L7|. The condition, beside the absence of infrared stable fixed points, is that the 
renormalization group trajectories cross the stability line, 


4(gi + # 2 ) + Pi(gi, 92 ) + @ 2 ( 91 , 92 ) — 0, 


a line at which the potential is minimized. Provided the potential vanishes at zero field, 
this should happen in the region where the potential is negative, gi + g 2 < 0. For N > 3 
the trajectories generated by the beta functions of (| 8 |) do cross the line of stability. The 
conclusion of Pisarski is that there is a first order phase transition for Nf > 3, and parity 
conserving mass is generated spontaneously for all Nf. 

This result seems incompatible with the picture of a chiral phase transition at some 
critical value of N c ~ Nf which we have described in the previous sections. In the remainder 
of this Section we shall list some of the possible reasons for this apparent lack of agreement. 


• In the conventional applications of the 4 — e-expansion to study the critical behavior 
of a finite temperature 4 dimensional theory, the fermions cannot contribute to critical 
fluctuations because they have antisymmetric boundary conditions in the time direction 
and thus their ‘mass’ is 0(1//?). Therefore, fermionic degrees of freedom need not be 
taken into account in the computation of the /3-function. In the zero temperature 
3 dimensional theories which we are studying here, however, the mass of fermions 
vanishes at the chiral transition point, thus their contribution is important and should 
also be taken into account. Thus if the theory has only four fermi interactions and no 
gauge interactions it belongs to the universality class of Yukawa theories. Therefore to 
do a reliable e expansion analysis one should add an appropriate Yukawa coupling to the 
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scalar theory in equation (7) and look for fixed points of the resulting beta functions. 
In the case of QCD3 the situation is somewhat complicated. If confinement effects 
are small , one must also include gauge interactions in the computation of the beta 
functions. It has been conjectured that when the number of fermion flavors exceeds 
a critical value .the long range attractive gauge interactions get screened leading not 
only to chiral symmetry restoration but also to deconfinement citeappelquist2. If this 
is the case Pisarski’s calculation should be appropriately modified to take into account 
the added complications. A calculation of effective scalar potential for QCD3 including 
gauge and fermion contributions has been carried out in the large flavor limit Jl^|. It 
is seen that the effective potential shows non analytic behavior as the critical point is 
approached from the symmetric phase. 


• The symmetries of the effective action (^) are not quite the same as that of three dimen¬ 
sional QCD described by (fj). Action m has a reflection symmetry, x which 

is distinct from parity and which is absent from the original gauge theory. This sym¬ 
metry cannot be broken by relevant terms in the framework of the Landau-Ginzburg 
approach. It is broken by the topological term in the sigma model (f|). However, it 
is impossible to take into account the influence of the topological term in the epsilon 
expansion since it cannot be continued in dimensions. 


There are known examples of non-perturbative infrared stable fixed points. Even in the 
scenario of the 4 — e-expansion starting with a Landau-Ginzburg effective action, that 
is valid in the large N c limit, it is very possible that there exists a non-perturbative 
infrared stable fixed point. Such fixed points are known to occur in 3-dimensional 
models with fermions and four-fermion interactions B. |19|. For example, consider pO 


S = J d 3 x (V(7 ' d + 0)^ + TT^ 2 ) ( 9 ) 

Eliminating <f> by using its equation of motion, produces the 4-fermi coupling 

which is invariant under a global SU(N) flavor symmetry. This model is renormaliz- 

able in the large N expansion. Integrating the Fermions gives the effective scalar field 

theory 


S'eff = —iVTrln( 7 - d + </>) + 

8 


N A 


2A 


2 


( 10 ) 













for which 1/iV is the coupling constant. It is known that this model has a non- 
perturbative fixed point at large A. Note the difference between this model and gauge 
theory. Since the former is strictly renormalizable, the ultraviolet cut-off A appears as 
a dimensionful parameter in addition to the dimensionful coupling A, thus providing 
an additional parameter AA. It is this parameter which is tuned to obtain critical 
behavior. 


The 4 —e expansion is known not to be reliable in some models. The Higgs-gauge model 
with N complex scalars, which is a model for superconductors, was found by Halperin, 
Lubensky, and Ma |2T|] to have a first order phase transition, due to the lack of infrared 
stable fixed points. This is a generalization of the Coleman-Weinberg argument. It 
is known that this system in the Type II regime (roughly speaking the Higgs mass is 
larger than the gauge mass) has the same critical properties as the smectic-nematic 
phase transition in liquid crystals. Experimentally, the latter system was found to 
have a second order transition with XY exponents. The problem has been studied by 


Monte Carlo simulations [22]. It was shown that at e = 1 there is a new fixed point 


with XY exponents, which is absent at e = 0. March-Russell [|23]| came to the same 
conclusion by studying the gauge interactions of N complex p-vector Higgs bosons in 
U (p)-gauge theory and comparing 4 — e and 2 + e expansions. In general, sigma models 
have been studied in the 2 + e expansion by Brezin, Hikami, and Zinn-Justin. They 
found fixed points and nonzero coupling for e > 0. This is in contradiction with results 
obtained when one investigates linear sigma models in 4 — e expansion, though the 
physical degrees of freedom driving the phase transition and the universality class are 
believed to coincide in the two models. 


If, as we have noted above is possible, the analysis of the 4 — e expansion is reliable in the 
limit of large N c , it is still implies that the chiral phase transition is a fluctuation induced 
first order transition. This comes about by the following reasoning: In the zero temperature 
gauge theory defined by action ([!]), the parameters are Nf, and N c . 

The coefficients of the operators in the effective Landau-Ginzburg theory (|7|) may depend 
on the parameters of the model. In the case of QCD in 2+1 dimensions, the only relevant 
parameters are N c and Nf. There is a phase transition if it of possible to tune N c and Nf so 
that the quadratic terms in the Landau-Ginsberg potential vanish. This would give a critical 
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curve in the N c — Nf- plane. The results of [h] imply that this curve is approximately linear, 
Nf ~ N c . The question of whether the transition is second order or a fluctuation induced first 
order transition remains. It is not clear that the analyses of the phase transition described 
in Section 3 distinguish a second order transition from one which is weakly first order. This 
could in principle be resolved using an improved 4 — e expansion. Also, the model which we 
construct in the next Section has the same flavor symmetries as 2+1-dimensional QCD and 
has a second order transition. If they indeed lie in the same universality class, this would 
indicate that the chiral transition at Nf ~ N c in QCD is indeed second order. 

5 A gauge model with global symmetry U(Nf) v 

In this Section we shall construct a model which has the same global symmetries as QCD, 
is renormalizable in the large N expansion, and has a second order chiral phase transition. 
Since the model shares the Landau-Ginzburg potential with QCD, claims of a proof of first 
order phase transition based on the Landau-Ginzburg potential and the 4 — e-expansion 
are invalidated. The model will consist of four-fermion interactions which drive the phase 
transition. To illustrate, we begin with a simpler version, a theory with only four fermion 
interactions and where the color symmetry is not gauged. 

5.1 A toy model with 4-Fermion interaction 

Let us begin this section with a toy model with 4-fermion interactions with Nf “flavors” and 
N c “colors” of fermions -0“ a = 1. .. N c , a = 1. .. Nf and an action which is a generalization 
of (D 

•S' = J d 3 x ($(7 ■ <9 + </>)V> + “7^ tr (^ 2 )) • ( n ) 

Here <p a p is a traceless Nf x Nf Hermitian matrix. Integrating over (j) gives the four-fermion 
interaction 

- J <^x^ a T A i\) a ^ h T A ^ b , ( 12 ) 

where both flavor and color indices have been omitted. T A is the generator of the flavor 
symmetry group, SU(Nf). One could modify interaction term (|l^) by including an isoscalar 
four-fermi interaction, (?/^) 2 . That case will be discussed briefly later. 

This model can be readily analyzed in the large N c expansion. If the color group were 
gauged, the necessity of summing planar diagrams would render this expansion intractable. 
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In the following section we shall suggest a way to overcome this difficulty. In this section we 
shall study the model (|TTD . 

The fermion fields can be integrated out to give in leading order of N c the following 
effective potential for the scalar held 


V[4>] = N c tr 






(13) 


where the trace is over “flavor” indices. This potential predicts a second order phase tran¬ 
sition at the point where the effective mass parameter 


t_ A A 


A 


TT* 


(14) 


vanishes. When t < 0 , the held 0 acquires a vacuum expectation value, breaking the U(Nf ) 
symmetry. The tree approximation gives the mean held value, 0, for the exponent 7 , dehned 
by 


R crit [0] ~tr|0| d(1 +7). 


(15) 


As we can see there is still no dependence of the universality class on N c at tree level. 

The effect of higher order corrections will be sketched below. In 0(1/N C ) one has to 
calculate contributions coming from the huctuations of the helds. The corrected potential 
has the form 


V[4>\ = N c 


tr0 2 -|-tr(0 2 ) 3 / 2 + —TRln A 


67T 


N r . 


where TR is a trace over spacetime as well as flavor indices and 


(16) 


\jki(x,y) = j \ 5ii5 jk - -^-SijSklJ S(x - y) + Tr(x| -. q + ^ \v)M ^, , g + ^ \ x )jk, (17) 

is the self-energy of the scalar to order IV” 1 . 

The global havor symmetry can be used to diagonalize the order parameter 0. Here, we 
shall assume the symmetry breaking pattern (when Nf is even) which preserves parity and 
an U(Nf/ 2) x U(Nf/ 2) havor symmetry 


= diagonal (m, —m, —tn ,..., — m). 


Then, the effective potential has the form 


N c Nf 9 iv c jvf o aiv f „ in 
V e s = —-— -tm 2 H- -rrf -|-f mr In —. 


N r N 


2 N3 


6 n 


37r 3 




(18) 


(19) 
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The cubic term of the effective potential comes from the fermion ring vacuum diagram, 
while the last term from the boson ring vacuum diagram. The effective mass parameter t 
also renormalizes at this order. We have also cancelled a logarithmic ultraviolet divergence 
by a wave-function renormalization for the scalar field, m —> zm and appropriate choice of 
z. 


The logarithmic term in the effective potential is reminiscent of a similar term in the 
Coleman-Weinberg [T^] effective potential in four dimensional theories. Furthermore, as in 
the Coleman-Weinberg analysis, the critical model (at t = 0) indeed has a minimum away 
from the origin and an apparent first order phase transition with the exponentially small 
order parameter 


7T 2 N r 


mo = yuexp 


4 N 


f, 


( 20 ) 


However, as in the Coleman-Weinberg scenario with one coupling constant, this minimum is 
at a point outside of the perturbative regime - it is for an exponentially small order parameter 
and the logarithm in the effective potential invalidates the perturbative result in the region 
where the order parameter is close to zero. 

As usual, the trouble with the effective potential near m = 0 can be fixed using the renor¬ 
malization group. This amounts to understanding how higher orders in logarithms would 
contribute to the effective potential. For the present case (unlike in the four-dimensional 
model considered by Coleman and Weinberg) we shall argue in the Appendix that the cor¬ 
rections exponentiate to produce an effective potential of the form 


Wff = 


NcNf tm 2 + NcNf 
2 6t r 


( 21 ) 


This improved effective potential exhibits a second order phase transition at t — 0. 

The exponent 7 defined in (pA|) is 0(1/ N c ) , 7 = ANf / (37r 2 N c ). The exponent is identical 
with the anomalous dimension of the boson. I 11 other words, the universality class depends 
on N c . In the broken phase, / 1 2 < 0, the eigenvalues of matrix (cp) are m, = ±mo, where 
777,0 ~ (Itl) 1 ^ 1 " 1 " 37 ^. 

It is worth mentioning what happens if a flavor neutral term is included in the four- 
fermion interaction. Then in addition to the adjoint representation scalar boson (isotensor 
boson), one is forced to introduce a flavor neutral (isoscalar) boson as well. The isotensor 
and isoscalar bosons have, in general, different couplings to fermion pairs. On one hand, 
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new divergent diagrams appear and both couplings start to run. The anomalous dimension 
of bosons is not any more twice the anomalous dimension of fermions. On the other hand, 
some crucial features of the pure isotensor theory survive. The logarithmic derivatives of 
couplings and the anomalous dimensions are independent of the couplings constants; they 
depend on Nf and N c only. The running fermion-boson coupling constant exponentiate into 
forms g ~ rri x ( N f)/ Nc . In a similar manner, the three boson couplings exponentiate as well. 
Thus, even in this case, the effective potential exponentiates into a pure Nf and N c dependent 
power of the the expectation value of the field, m. The phase transition is of second order 
and critical exponents are dependent on Nf and on N c . Details of the discussion of this 
model will be given elsewhere. 

5.2 Gauging the “color” group: Chern-Simons gluo-dynamics 

The four-fcrmi theory which we have been discussing has a second order phase transition. 
However, the flavor symmetries is very different from that of 2+1-dimensional QCD. One 
can come closer to QCD and remove the color degeneracy of the model by gauging the color 
symmetry. We shall introduce gauge fields in such a way as to preserve, as much as possible, 
the solvability of the four-fermion model in the large N c limit. 

Consider the action 

S — J 'Px +7 ■ D a + 4> + + x(l ■ Db + <P ~ M)x 

( 22 ) 

where Da = d + iA and Db = d + iB. Here, we have assumed that there are an even number 
of colors and have gauged the group SU(N c /2 ) x SU(N c /2). As a kinetic term for the gauge 
field, we have introduced the Chern-Simons terms with a dimensionless coupling constant 
g. This makes the gluon non-propagating. It is very important to observe that even though 
there are Chern-Simons terms in the theory, it is still parity invariant. In fact, this theory 
has the same flavor symmetries as QCD: it is parity symmetric, the parity transformation, 
besides the obvious spacetime transformation, also interchanges the two types of fermions x 
and ^ and the two types of gauge fields A and B. This symmetry insures that if there is no 
bare fermion mass it will not be generated to any finite order in perturbation theory. 

This theory also has the flavor symmetry U(Nf) and one can study the possibility of a 


N C A 
' 2A 


tr(0 2 ) + 


Nc 

2 g 2 


-tr 


Ad A + -A 3 - BdB 
3 
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phase transition with symmetry breaking pattern U(Nf ) — » U(Nf/2) x U(Nf/2). Its one 
large difference with QCD is that here the fermions have a bare mass, M. They are not 
massless at the critical point unless we also tune M — 0. This model thus allows us to 
explore two possibilities, one which is like QCD with M — 0 and where massless fermions 
are important to the critical behavior and the other where M/ 0, fermions are all massive 
in the region of 0 ~ 0 and one would expect that perhaps the Landau-Ginsburg theory and 
epsilon expansion gave an accurate description of the fluctuating degrees of freedom. 

For the gauge interactions, to obtain the large N c limit one must sum all of the planar 
diagrams which contribute to the effective potential. The leading order is N 2 and contains no 
fermion loops. It therefore does not couple to the order parameter 0 and does not contribute 
to the effective potential for 0. The next to leading order, ~ N c , contains one fermion 
loop, no scalar lines and all insertions of planar gluon lines. The problem of summing all 
contributions of planar gluons is as yet intractable. Here, we shall remedy this situation 
by using a simultaneous large N c and small g 2 expansion p5[. In this expansion, unlike 
conventional gluo-dynamics where the gluons would have a Yang-Mills action, the Chern- 
Simons gluo-dynamics is infrared finite, i.e. infrared divergences do not ruin the perturbative 
expansion 0- 

Then, when M — 0 (the QCD case), the order g 2 N® contribution to the effective potential 
vanishes. Then, in this model, the effective potential is given by (jT^) with corrections being 
of order N^g 2p with k + p > 1. Furthermore, since the gauge interactions do not contribute 
to the logarithmic divergence to the order we have considered, the renormalization group 


analysis as well as the result (|2l|) are also valid. Thus, this theory has a second order 
phase transition and has the same flavor symmetry as QCD. 

The analysis of this problem when M / 0 is more difficult and thus far is done only to 
the leading order. There, the effective potential has the form 


VI 


Mr 
eff [ 


= N c tr 


W 1 + w + (! + £)£ 

|M 2 + ^|0| + i|0| 3 


|0| < M 
101 > M 


(23) 


This system has a marginally first order behavior. When t > t c — —Mj2ir the minimum is 
at 0 = 0. Then at t — t c the effective potential has zero curvature in the region |0| < M 
and is convex where |0| > M. When t > t c the minimum begins at |0| = M and increases 
as t decreases. This is very similar to a first order behavior: the order parameter 0 jumps 
from 0 = 0 to |0| = M at t = t c . 
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It is intriguing that this “almost first order” behavior occurring for this case is probably 
closer to the prediction of the epsilon expansion than the case of pure QCD. However, even 
here we conjecture that higher order corrections smooth the behavior of the effective potential 
and the phase transition remains second order. Detailed behavior of this model in the next 
order of expansion are interesting and are the subject of further study. 

6 Conclusion 

In this Paper, we have analyzed a pure four-fermion model and a gauged four-fermion model 
in the large Nq limit. We have shown that both of these models exhibit second order phase 
transitions. 

Perturbatively, the models appear to have fluctuation induced first order transitions. 
However, we have shown that higher order corrections lead to exponentiation of logarithmic 
terms and restoration of second order behavior. 

In the gauged four-fcrmi model, we could consider two cases - one where the fermions 
are massless at the phase transition and the other where they are massive. In the case 
where they are massless, at least to the order of perturbation theory which we considered, 
the nature of the phase transition was identical to that in the pure four-fermion theory. In 
the case where the fermions are massive, however, the behavior was markedly different - it 
exhibited an “almost first order” phase transition. It is intriguing that the 4 — e expansion 
predicts a fluctuation induced first order transition for this model. It is also interesting that 
as the fermion bare mass goes to zero, the transition becomes second order (the jump in the 
order paramater as the coupling is varied through the critical point is proportional to the 
bare mass). 


A Appendix 


The exponentiation of the potential is proven by the method of renormalization group. The 
held theory we consider is defined by action (|TT|). The fermion bubble corrections are summed 
to give a modified boson propagator, which has the form 


D~ l 


*? 


(P) = 


A A \ m p 2 + {rrii + rrij) 2 _ 1 p 
A 7T 2 / 2tt + Anp 2 m 
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1 

2 


+ 0(m 4 p 3 ),(24) 


/A A \ p (m ? ; + 77lj) 2 

V A 7r 2 / 8 8p 


m[3(mj + rrij) 2 — 4m 2 ] 
bnp 2 


where m* = ±m is the mass of the ith fermion. 

The renormalization group functions can easily be calculated in 0(Nf/N c ). The only 
primitive divergent diagram is the fermion self-energy correction diagram in this order. Thus, 
the coupling constant does not run. The logarithmically divergent part of the fermion self 
energy diagram is 

4 Mf A 

A S- F \k)^—f- 2 (25) 

3N c n z /j 

The boson propagator gets logarithmic corrections only due to the dissociation of the 
boson into two fermions. Thus, the leading order anomalous dimension of the boson is twice 
the anomalous dimension of the fermion. In other words, as shown by (f25P, 


lb 


4 N f 
3N c n 2 ' 


(26) 


At the critical point, the Callan-Symanzik equation 


' d d \ 

s +7tm *s) 


V{m) = 0 


(27) 


has the boundary condition 


lim 


VA 


= - 771 . 


N f /N c ^o N c Nf 6n 

Here t = log/i/m. Then the solution of the Callan-Symanzik equation is straightforward. 


(28) 


Ktr = — m? exp 
bn 


fa' 

/o 1 - 71.(7') / 


(29) 


Substituting (^) into (|29|) gives exponentiated form (pT|) . 
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